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THERAPEOTIC DECISIONS SYSTEMS AND METHOD USING STOCHASTIC 

TECHNIQUES 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] This application claims benefit of U.S. 
Provisional Application No. 60/299,578, filed June 19, 
2001. 

BACKGROUND 

[0002] The determination of how to treat a highly ill 
patient can be very difficult, especially under the 
chaotic and emotional conditions that may exist in the 
conditions surrounding around such a highly ill patient. 
For example, such an ill patient may be treated in the 
context of an emergency situation such as shock, high 
risk surgery, trauma and other acute conditions requiring 
emergency treatment. It has been found by the inventors 
that diagnosis errors commonly occur under those 
emergency conditions. 

[0003] It is also difficult, especially in these 
emergency situations, to evaluate the timing of a given 
therapy. A given therapy that improves the outcome at 
one point may actual produce harm at a different point, 
or when used in the wrong amount for the wrong situation. 
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For example, fluid therapy used at the initial phase of 
resuscitation may be critical. However, use of excess 
fluids may lead to pulmonary edema and cardiac failure. 
Such errors can be lethal in emergency situations. 
Moreover, many different injuries, shock being one, may 
be easy to recognize in late stage when therapy is often 
ineffective. If hypovolemia is inferred from 
tachycardia, hypotension, and falling hematocrit, these 
superficial manifestations of shock may be overtly 
corrected by transfusions, fluids, and vasopressors, but 
still without adequately restoring the underlying 
circulatory functions. 

[0004] It has been found that effective resuscitation 
of acute life-threatening emergencies achieves optimal 
physiological goals as early as possible. When 
circulatory mechanisms are identified earlier and treated 
more vigorously to specified physiologic target goals, 
outcomes are improved. 

[0005] The incidence, mortality, morbidity, and costs 
of life threatening illnesses and injuries are 
extraordinary. Among the 1.9 million deaths annually, 
about half are from acute illness associated with shock 
and lethal organ failure. There are over 34 million 
surgical operations annually in the U.S., with an. overall 
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mortality of 1%, but high-risk surgical patients have 
mortality between 25 and 33%. Postoperative deaths are 
often due to adult respiratory distress syndrome (ARDS) , 
which has an incidence of 150,000, a mortality of around 
40%, consumes an average of two weeks in the ICU, and in 
this country costs about $1,950,000,000 annually. There 
are about one-half million septic patients, of whom 40% 
develop shock, with 50 % mortality. Septic shock is the 
13th leading cause of death, and the most common cause of 
ICU deaths. Hospital costs for septic shock are over 5 
billion dollars annually. 

[0006] Conventionally, shock is classified as 
hemorrhagic, traumatic, postoperative, neurogenic, and 
distributive or septic. Usually this classification is 
applied to the later, fully developed, clinical syndromes 
of shock. Analysis of circulatory mechanisms when 
patients are admitted to the ICU in the late stage, after 
organ failure has occurred, is extremely complex because 
of the many interacting clinical, physiological, and 
immunochemical problems. Over 12 separate cascades of 
chemical mechanisms have been described. Many of these . 
have been considered to be the cause of shock, and major 
efforts have been expended to reverse them. Of the 11 or 
more large scale, multicenter, randomized clinical 
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trials, only the recombinant human activated protein C 
(drotrecogin alfa activated) has shown inqprovement in 
mortality, which decreased from 30.8% in the placebo 
control group to 24.7% in the protocol group, a decrease 
of 6.1% (38). The many interacting immunochemical 
mechanisms make the problems of sorting out causal 
relationships more difficult, and the therapy less 
effective, after organ failure and sepsis have become 
established in the late stages of critical illness. 
[0007] There are six commonly used outcome predictors: 
the Acute Physiology and Chronic Health Evaluation 
(APACHE), Glasgow Coma Scale (GCS) , the penetrating 
abdominal injury (PATI) score, the Therapeutic 
Intervention Scoring System (TISS) , The Revised Trauma 
Score and Injury Severity Score (TRISS) and the Trauma 
score. All of them assess categories of patients in terms 
of probable mortality. The are not intended to predict 
mortality risk for a specific single patient. Those 
systems are rarely used for guiding actual patient care. 
They find their grSeatest use as a tool for administrative 
and management research and decisions. Longitudinal 
scoring on successive days may reflect continuous 
improvement, but has no physiologic or therapeutic 
relevance for any specific individual patient. Moreover, 
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none of them identify specific underlying physiologic 
mechanisms or problems as does the present approach. In 
addition, none of them considers specific therapeutic 
recommendations to be suggested, real-time, to the 
physician. Finally, none are able to recommend titration 
of therapy to alleviate the underlying hemodynamic 
problem or to achieve optimal goals. 

BRIEF DESCRIPTION OF THE DEU^WINGS 
[0008] These and other aspects will now be described 
in detail with reference to the accompanying drawings, 
wherein: 

[0009] Figure 1 shows a block diagram of a patient 
being monitored, and the system described herein; 
[0010] Figure 2 shows a flowchart of the basic 
operation of the system; and 

[0011] Figure 3 shows a flowchart of relationship with 
the database. 

SUMMARY 

[0012] The present technique, in contrast to the 
above, defines techniques for detecting, and optimally 
correcting, for each individual, patient, in a real-time 
manner, the initiating hemodynamic mechanisms of shock 
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and circulatory dysfunction, in conjunction with a new 
outcome predictor and real-^time guidance and decision 
support for optimizing therapy. This is very different 
from the above scoring systems currently in use, which 
provide only reports concerning prognosis. 
[0013] The present invention teaches a system and 
machine of patient monitoring of various parameters and 
decision support that is guided by a stochastically 
analyzed operation. This operation can evaluate and 
optimize the therapy that is used especially in such 
difficult situations. 

[0014] Another aspect is an outcome predictor that 
calculates the probability of survival for any individual 
patient. This predictor, in an embodiment, is based on a 
program of stochastic/probability analysis developed from 
a database that includes clinical conditions. 
[0015] An embodiment determines the patients state (S) 
using a specified mathematical criteria, and determines 
patients in the database most similar to the patients 
state. Individualized therapy that maximizes the 
probability of survival may be determined from those 
other patient states that are closest to the patient' s 
state. 
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DESCRIPTION 

[0016] Many or most deaths in acute illnesses are 
based on circulatory factors. An aspect of this 
application describes characterizing a patient based on 
hemodynamic parameters and comparing those parameters to 
previously obtained information in a database. While the 
present embodiment describes hemodynamic parameters^ it 
should be understood that other parameters besides 
hemodynamic parameters can be analyzed. Hemodynamic 
measurements, however, may be specifically advantageous 
since many previous studies and clinical trials have 
established relationships between such hemodynamic 
parameters and patient results (where patient results may 
include survival, long term injury, or other results) . 
Also, noninvasive monitoring technology is available in 
many different forms for obtaining hemodynamic 
measurements. 

[0017] An embodiment shown in figure 1. A monitoring 
part 100 monitors various parameters of the patient. The 
monitoring may be both noninvasive and invasive, although 
noninvasive monitoring may be preferred. As described 
herein, the monitor may monitors various characteristics 
including cardiac index, blood pressure, arterial and 
cutaneous oxygen, carbon dioxide level, and other 
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parameters as discussed. Multiple different parameters 
may be monitored. The output parameters 110 represents 
all of these different parameters. Moreover, while the 
monitoring element 100 is shown as being in a single 
location, it should be understood that the monitoring may 
be carried out at multiple different locations. 
[0018] The output parameters 110 are coupled to a 
processor 120, which receives the data and carries out 
operations using the data as described herein. The 
processor 120 is also associated with a database 125. 
The database includes various information about the 
patients that have been treated in the past. Again, 
further details on this are described herein. Both the 
processor 120 and database 125 can be either local, or 
may be remote from the patient, in which case the link 
that carries the data 110 may be a network connection 
e.g. the Internet. 

[0019] The database includes information about other 
patients as a function of all or some of the monitored 
variables. The database and processor determines 
mathematically the state of this particular patient by 
finding the mathematical nearest neighbors. These 
mathematical nearest neighbors can also be analyzed to 
determine what percentage .of these patients live or die. 
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and/or how successful specified treatments on these kinds 
of patients have been in the past. This can be used, in 
turn, to make decisions about what kind of therapy might 
be the best kind of therapy both in the short-term and in 
the long term. 

[0020] All measurements are preferably started with 
the patient's admission to the emergency department (ED), 
When the paramedics radio the ED that a trauma patient is 
coming in, the trauma team on call is alerted and is 
ready in the ED to receive the patient. The patient is 
then followed to the radiology department, the operating 
room (OR), the intensive care unit (ICU) , and ultimately 
to the patient's survival or death at hospital discharge. 
In surgical patients, the noninvasive hemodynamic 
measurements are started with onset of anesthesia, 
continued throughout surgery, and followed in the ICU, 
step-down units, and hospital floors. 
[0021] Monitoring element 100 may monitor the 
following: 

[0022] -Bio impedance. This is a device that measures 
the biological impedance of the user. Various injuries, 
such as chest wall edema or chest tubes may cause the 
electrical signals to be distorted, but compensation for 
this distortion can be made. 
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[0023] Biological impedance may be measured using an 
improved continuous noninvasive impedance system (IQ 
device) has developed by Wang et al at Drexel University 
and marketed by Wantagh, Inc., Bristol PA. They use 
noninvasive disposable prewired hydrogen electrodes 
positioned on the skin, three EKG leads placed across the 
precordium and left shoulder, and a 100 kHz, 4 mA 
alternating current passed through the patient's thorax 
by the outer pairs of electrodes. The low voltage is 
sensed by the inner pairs of electrodes. The voltage 
sensed by the inner electrodes captures the baseline 
impedance (Zo) , the first derivative of the impedance 
waveform (dZ/dt) , and the EKG. The EKG and the 
bioimpedance signals are filtered with an all-integer- 
coefficient method that simplifies computations and 
decreases signal processing time. 

[0024] The signal processing technique is based on 
time-frequency distribution technology that provides high 
signal-to-noise ratios that increases the speed of signal 
processing and permits essentially real time operation 
without deterioration of accuracy. It shows the 
relationship between the time, frequency, and power of 
the signal that would otherwise be unobtainable from the 
standard fast Fourier transform. This is particularly 
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important in identification of opening and closing of the 
aortic valve ("B" and "C" points), which are often buried 
in artifacts that make them difficult to identify. 
Variables available for display and printout include: 
cardiac output and index, stroke volume, Zo (baseline 
impedance), MAP (mean arterial pressure), HR (heart 
rate), systemic vascular resistance, dZ/dt max (maximum 
increase in impedance rate), pre-ejection period, 
ventricular ejection time, acceleration index, dz2/dt2, 
and end-diastolic volume. 

[0025] -Pulse oximetry (machine from Nellcor, 
Pleasanton, CA) continuously measures arterial oxygen 
saturation (Sa02) . Appreciable changes will be noted and 
confirmed by the standard in vitro blood gas analysis. 
Pulse oximetry has been extensively evaluated and is 
routinely used. The pulse oximetry may be used to 
[0026] Compare Pulse Oximetry with In Vitro Hemoglobin 
Oxygen Saturation, and to 

[0027] Compare Transcutaneous Oxygen Tension (PtcOa) 
and Oxygen Consumption. Oxygen consumption (VO2) values 
may be compared with simultaneous Ptc02 at the initial 
baseline period, the nadir, and after resuscitation. The 
Ptc02 nadir occurred an average of 12.1 + 7.7 min (p<.05) 
before the VO2 nadir, indicating that while both methods 
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reflect tissue perfusion/ PtcOa decreases occur firsts and 
provide an earlier warning of inadequate tissue 
perfusion. 

[0028] -Survivor and Nonsurvivor Patterns in High-Risk 
Surgery and Trauma. 

[0029] -Therapeutic Effects of Blood, Albumin, 
Crystalloids, Dopamine, and Dobutamine, evaluated by 
Invasive Monitoring. 

[0030] -Noninvasive Monitoring of the Effects of 
Packed Red Cells, Albumin, and Crystalloids. 
[0031] -Comparison of Temporal Patterns of Hemodynamic 
Data with Immunochemical Mediators. 

[0032] -Transcutaneous Oxygen and Carbon Dioxide for 
Estimation of Tissue Perfusion. 

[0033] Standard transcutaneous oxygen tension {Ptc02) 
measurements (Novametrix Medical Systems, Inc., 
Wallingford CT) will be continuously monitored throughout 
the observation period. Previous studies have shown well 
the ability of Ptc02 to reflect tissue oxygen tension. 
Ptc02 reflects delivery of oxygen to the local area of 
skin; it also parallels the mixed venous oxygen tension 
except under late or terminal conditions where peripheral 
shunting may lead to high mixed venous hemoglobin 
saturation (Sv02) values (18). It uses the same Clark 
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polarographic oxygen electrode routinely used in standard 
blood gas analyses (15,18-20). Oxygen tensions are 
determined in a representative area of the skin surface 
heated to 440C to increase emissivity of oxygen across 
the stratum corneum, and to prevent vasoconstriction in 
the local area being measured (19) • Transcutaneous C02 
monitoring (Novametrix Medical Systems, Inc., Wallingford 
CT) of the skin surface will be continuously monitored by 
the standard Stowe-Severinghaus electrode (14) in the 
same sensor unit. 

[0034] Covariate information may also be obtained. 
[0035] The following 30 covariates may be measured: 1) 
age, 2) gender, 3) preoperative, intra-operative, and 
postoperative estimated blood loss, 4) blunt or 5) 
penetrating truncal trauma, 6) nontruncal (extremity) 
injury, 7) spinal cord injury, 8) blunt or 9) penetrating 
cardiac injury, 10) pulmonary contusion, 11) pelvic 
fracture, 12) long bone fractures, 13) head injury, 14) 
brain death, 15) early stage (<12 h) , 16) middle stage 
(12 to 24 h), 17) late (>24 h) stage, 18) cardiac 
insufficiency (reduced cardiac reserve capacity 
determined by responses to standardized doses of 
transfusions, and fluid challenges), 19) bacterial 
contamination, sepsis or systemic immune response system 
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(SIRS), 20) respiratory dysfunction or failure 
immediately prior to the present acute illness, 21) pre- 
illness renal insufficiency or failure, 22) pre-illness 
hepatic failure, 23) nutritional insufficiency or 
failure, 24) uncontrolled diabetes, 25) pre-illness 
essential hypertension, 26) cardiac injury, blunt or 
penetrating, 27) cardiac arrest, 28) pregnancy, 29) 
Glasgow coma score, and 30) the injury severity score 
(ISS) . 

[0036] Quantitative Assessment of Continuously 
Monitored Noninvasive Variables as Net Cumulative Amount 
of Excess or Deficit may be an important factor- While 
monitored values can be directly observed, many show 
considerable variability, which may obscure the 
underlying pattern. To overcome this problem, we 
calculate the net cumulative excess or deficit of 
physiological variables by integrating the areas between 
the curves of monitored data and selected normal or 
"optimal" values. This provides a quantitative measure of 
the overall deficit in cardiac, pulmonary, and tissue 
perfusion function. This may provide a more useful early 
evaluation of the patient's physiological responses than 
simply following empirical scoring systems (Apache, 
etc.,) which are designed only for prognostic reports. 
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Real-time suggestions concerning the best therapeutic 
course of action for the patient are simply beyond their 
scope. 

[0037] The processor 120 may carry out a program that 
is described in the flowchart of figure 2. At 200, all 
of the variables described above may be obtained, as well 
as other variables which are described herein. These 
variables are preferably obtained by noninvasive 
techniques, although they can alternatively be obtained 
by more invasive techniques. At 205, each of these 
variables is evaluated. The evaluation may be carried 
^out by comparing the variable with its normal or optimum 
value and a temporal pattern during the observation. The 
variable may be mathematically integrated over time 
versus normal or optimum variables. Net cumulative 
deficits or excesses are then calculated for each 
individual patient. 

[0038] Table 1 shows the different variables and how 
they compared in some survivors and non survivors. The 
variables obtained include cardiac index, mean arterial 
pressure, arterial hemoglobin saturation and 
transcutaneous oxygen. 
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[0039] The mean values similarly show this trend: 
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[0040] In an experiment, out of 103 patients who 
survived their illnesses, 98 were correctly predicted to 
survive the illness by comparing their values with other 
similar patients- Of the 48 who did not survive, 30 were 
correctly predicted not to survive. This compares with 
previous studies in which much poorer results were 
obtained. Thus far, using this approach those with 
outcome probability is greater than 80 percent have 
survived, and always outcome probability is less than 35 
percent have died. The current in. 

[0041] In general, the specific monitored goals of 
therapy will be based on the current empirical findings 
of the survivors present in our current patient database: 
cardiac index > 4.5 L/min/m^; systolic blood pressure > 
120 mmHg; pulse oximetry > 96%; transcutaneous oxygen 
Ptc02/Fi02 ratio >200; heart rate < 100 beats/min. Packed 
red cell transfusions will be given to maintain Hct > 30% 
or hemoglobin concentrations >10 g/dL. The time to 
obtaining these goals, the amount and type of fluids will 
be noted. 

[0042] At 210, an outcome predictor is used with the 
variables. The outcome predictor uses the probability of 
survival as being a measure of the severity of illness 
before and after treatment of hemodynamic problems that 
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may contribute to the non survival pattern. Correction 
or failure to correct each mechanism by the specific 
therapy at each temporal stage during the acute illness 
is taken as a measure of the reversibility of the injury 
mechanism. Then^ the predicted changes in outcome are 
compared with actual survival or death at hospital 
discharge. The specific flowchart of the outcome 
predictor is shown in figure 3. The flowchart operates on 
dynamic program based on the "state" of the patient at 
any time, shown as the variable X, . Based on this state, 
an optimal therapeutic policy is determined that has, on 
the average, proven to work best for similar patients in 
similar conditions that are recorded in the database 125. 
The term * similar' here means a group of patients who 
have the same diagnoses who share the same set of 
specified covariates, and have similar hemodynamic 
states. 

[0043] This state of the patient at any time is 
therefore defined in terms of primary diagnoses, 
covariates, and hemodynamic measurements. x(t) is 
referred to as the state vector at time t. This is 
defined in terms of the hemodynamic measurements, their 
derivatives, and their integrals. 
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[0044] Assuming that there are L different kinds of 
measurements taken on a given patient. The set may 
include for example cardiac index, blood pressure, pulse 
oximetry, transcutaneous oxygen and CO2 tensions. Each 
measurement I is denoted as yL. The state vector is 
defined as a concatenation of the value yi, as well as its 
derivative or or each measurement type, denoted as y^f 
define the state vector as a concatenation of the value 
itself, its first and second derivatives y^, y\ , and 

its first integral \yidt ^ as follows: 



yx{t^yx{h\y\{t^^^ 



(2.1) 

i.e., for L different measurement types there will be 4L 
states. In practice, the derivatives and integrals are 
approximated by finite differences and siams of the time- 
ordered data of the database- Specifically, we will 
calculate the approximations. 
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]y,dt s y, {ft Xtk - tk-i) + (y,dt (2.2) 

0 0 

[0045] Since not all measurements are taken at all 
times in the database, smoothing or fitting the 
measurements by interpolation will be used to place them 
on a single common time grid, e.g., linear spline, cubic 
splines, etc. The best method for smoothing, beginning 
with the simplest linear splines, will be examined. 
[0046] The control input is one of M different control 
input that can be applied to the system. 

[0047] Specifically, the control input, m(/j) at time 7^ 
is assumed to be drawn from the finite set, 

"('*)€{« ,1/^} (2.3) 

[0048] The simultaneous administration of two or more 
types of therapy will be treated mathematically as a 
separate distinct mode of combination therapy. Hence M 
generally includes all single therapies plus relevant 
combinations of therapies, such as packed red cells and 
crystalloid fluids, fresh frozen plasma and transfusions. 
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[0049] It is convenient to think of the propagation of 
the patient's state at time , to his state x^^jat time 
tf^^^as obeying the following nonlinear dynamical system 
with process noise and parameters p, i.e.: 

4.i^f(4^^nP^^k) (2.4) 
[0050] For simplicity, p is discrete, and is assumed 
to be drawn from a finite set formed by enumerating all 
useful combinations of covariates, 

P^{pt>"^PM} (2.5) 
[0051] Both covariates and process noise help to 
explain the variability of patient responses seen in the 
database. The covariates help to distinguish gross 
differences in responses due to patients with major 
differences in the nature of their disorders and 
complications. Process noise helps to explain small 
differences between patients with the same covariates but 
different responses to the same therapy. It is a measure 
of unmodeled dynamics, or intra-individual variability, 
due to other sources of variability in the system. 
[0052] The outcome predictor calculates the value PC 

or probability of survival. 

[0053] For a given state x, the survival probability 
is denoted by S(k), which is calculated simply by 
extracting the N nearest neighbor states to x of patients 
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having the saine collection of covariates p (denoted as 
the set (xjj^l ^ N (x,p)) and by calculating the 
fraction of them that survived. For example, if W s of 
the N nearest neighbor states to x survived, then the PC 
is given as, 

S(x) - N s/N (2.6) 
[0054] This quantity is equivalent to the estimated 
survival probability, and is a useful measure of the 
severity of the patient's current illness. 
[0055] A stochastic controller is shown in figure 3, 
which does not require a model, but which rather uses the 
patient database directly . When a new patient is 
admitted, the analysis and control (decision support) is 
calculated by the following operations. 
[0056] Step 0. Take measurements at 3 distinct 
successive times to calculate a nominal state xj according 
to (2.2), and define a relevant covariate vector 
Define k=0 to initialize the algorithm at time zero. 
[0057] Step 1. Given the nominal state and a 

covariate vector /?° , extract the N nearest-neighbor 

states (denoted as {xj^^, = ^^^(x^,/?^)) that are closest to x^^in 

the database and share the same covariates. Here, a 
measure of ""closeness" is conveniently defined in terms 
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of the quadratic distance, d{x,xl)s{xl -xfw{xlx) 
(3.1) 

where the weighting matrix W is taken to be diagonal ^ 

W = diag[w^, (3.2) 

The best choice of weighting matrix W remains to be 
investigated further. 

[0058] Step 2. From the set of nearest-neighbors 
{jC;t}^_j calculate the fraction of them which survived, to 

define s(xl). 

[0059] Step 3. Sort the N nearest-neighbors states 
into m groups depending on the type of control input u, 
which was applied to it (as recorded in the database). 
This gives the group of sets, 

{uLu]l,J = l,..,.,m (3.3) 

[0060] For example, the subgroup {r];,w^}^^, denotes the 

subgroup of nearest-neighbor states which received a 
specific control (therapy) input of type i=l. Note that 

f^m^^N (3.4) 

since the state has a specific control associated with 
it. 
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[0061] Step 4. For each state in the i'th subgroup, 
calculate the subsequent state xf^, 

[0062] By moving to the next line in the database, 
i-e, xf^,=/(xif,i/„/7>J (3.5) 

Calculate the survival probability 5'(x^+,) associated with 
this new state (remember this requires extracting nearest 
neighbors to x^^j ) . Define the PC associated with the i'th 
control as 

A(".)=-£4L) (3.6) 

[0063] Step 5. Apply u to the patient, take new 
measurements (to define new state x^^^, ) , set A<-ife+l, and 
go to step 2 . ' 

[0064] An important feature of the above is the 
ability to defined the patient's unique "state" at any 
given time. I the database that is used preferably is a 
large database with multiple timelines describing 
clinical and hemodynamic patterns of subsets at different 
conditions. The states should be described in time that 
elapses from the onset of illness or from admission. 
[0065] Weighting criteria to specify the nearest 
neighbors have been based upon the clinical diagnosis, 
,co-morbid conditions, and hemodynamic patterns of 
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survivors vs. nonsurvivors. For example, for the series 
as a whole, discriminant analysis revealed significant 
differences between survivors' and ndnsurvivors' net 
cumulative deficits in Ptc02/FiO2 (p< 0.0001), SapOa 
(p<.005), CI (p<0.03), and the initial GCS value (p< 
.0001). Further, we will use, as weighting criteria, . the 
p values of differences in the temporal patterns of 
survivors' and nonsurvivors' values in each subset, for 
example, for subsets of patients with and without head 
injuries, blunt vs. penetrating trauma, truncal and 
nontruncal trauma, age stratifications, prior cardiac, 
respiratory, hepatic, and renal dysfunction or organ 
failure, etc. Finally, the hemodynamic responses to 
standardized test doses of whole blood or packed red cell 
transfusions given over a specified period of time 
(usually 1-hour) provide quantitative measures of each 
patient's cardiac reserve capacity. This is expressed as 
changes in cardiac index relative to corresponding change 
in PA occlusion (wedge) pressure or central venous 
pressure (Starling's myocardial performance curve) before 
and after a blood transfusion given over a one-hour 
period. Criteria for satisfactory vs. limited cardiac 
functional reserve capacity have been developed for each 
clinical subset and each time period. 
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[0066] In operation, the program has been used to 
evaluate individual patients states in terms of 
probability of survival based on initial baseline 
measurements. So far, all those with outcome 
probabilities > 80 percent have survived and all those 
with probabilities less than 30 percent have died. At 
300, the patients individual data is put into the 
database to determine theN nearest neighbors of the user. 
The survival probability is determined at 305, using the 
operator S, the patients are also sorted by control 
intervention at 310. The state response is propagated at 
315, and at 320, the survival probability is calculated 
conditioned on intervention. Stochastic probability 
control is determined at 325. 

[0067] The stochastic analysis and control program can 
be used for data mining to detect significant hemodynamic 
events and their relationship to therapy given. It will 
begin by tracing individual patient histories in the 
database through time, where at each instant, the ''state" 
of the patient (see definition earlier) is computed, as 
well as the patient's associated probability of survival 
(based on nearest neighbors) before and after specific 
therapy. A large change, either up or down, in survival 
probability indicates a potentially significant role of 
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the associated therapy. This flags the corresponding 
patient and time for a closer look at what is actually 
going on. The entire database will be processed this way^ 
and statistics about the various significant events will 
be compiled. Therapeutic interventions with the most 
efficacy or harm can be flagged for more detailed 
analysis of the contribution of each therapeutic modality 
and each hemodynamic variable to each patient's survival 
or death- . 

[0068] Although only a few embodiment have been 
disclosed in detail above, other modifications are 
possible. 



28 



wo 02/103492 



PCT/US02/19740 



WHAT IS CLAIMED IS: 

1. A system comprising: 

a monitoring device monitoring a plurality of 
physiological signals of a current patient; 

a database, which stores records indicative of a 
plurality of patients including at least some of said 
physiological signals for said plurality of patients and 
prognosis information for said plurality of patients; 

a processor, comparing said physiological signals to 
said stored records in a database, to determine 
mathematically patients in the database which are most 
similar to a current patient; and to make available 
information about treatment and results for said patients 
who are most similar to the current patient. 

2. A system as in claim 1, wherein said 
physiological signals include hemodynamic information. 

3. A system as in claim 2, wherein said hemodynamic 
information includes bio impedance information. 

4. A system as in claim 2, wherein said hemodynamic 
information includes temporal patterns within the 
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hemodynamic information. 

5. A system as in claim 2, further comprising 
obtaining covariate information indicative of a patient's 
individual characteristics. 

6. A system as in claim 5, wherein said processor 
mathematically determines said similar patients by 
determining a state variable for a patient based on at 
least said hemodynamic information and said covariate 
information. 

7. A system as in claim 6, wherein said state 
variable is defined as a vector over time, where each 
point in time includes measurements at that point in 
time. 

8. A system as in claim 1, wherein said state 
variable also includes derivatives of the hemodynamic 
measurements and integrals of the hemodynamic 
measurements. 

9. A system as in claim 2, wherein said hemodynamic 
measurements include at least cardiac index, blood 
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pressure, pulse oximetry information, and transcutaneous 
gas tension. 



10. A system as in claim 2, wherein said processor 
also carries out interpolation to place a plurality of 
said values on a common time grid. 

11. A system as in claim 1, wherein said state 
variable is also a function of non covariate parameters, 
and inconsistencies between said state variable are 
explained in terms of said covariate parameters. 

12- A system as in claim 7, wherein said processor 
models said patient as a control system, with said state 
parameter representing a patient's state, a control input 
representing a therapy, and the patient's state following 
a nonlinear dynamic system with process noise being 
explained in terms of covariate parameters. 

13. A system as in claim 12, wherein said processor 
models said patient with a survival probability based on 
nearest neighbor states to other patients in the 
database. 
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14. A system as in claim 7 wherein said processor 
uses said state variable to find nearest neighbors in the 
database which represent patients that are most similar 
to the current patient. 

15. A system as in claim 14, wherein said finding 
nearest neighbors comprises determining a quadratic 
distance between the current patient and other patients 
in the database. 

16. A system as in claim 14, wherein said processor 
uses a weighting matrix to determine said closest 
neighbors. 

17. A system as in claim 16, wherein said weighting 
matrix is based on relationships between hemodynamic 
patterns and survival versus non survival. 

18. A method, comprising: 

monitoring a plurality of physiological signals of a 

current patients- 
comparing said physiological signals of the current 

patient to stored physiological signals for previous 

patients stored in a database; 
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determining patients in said database who are most 

similar to the current patient; and 

determining information about treatments and results 

for said patients who are most similar to the current 

patient. 

19. A method as in claim 18, wherein, said 
determining information comprises determining a 
likelihood that the patient will survive. 

20. A method as in claim 18, wherein said 
determining information comprises determining a therapy 
for the current patient which has proved most 
statistically successful for said patients who are most 
similar to the current patient. 

21- A system as in claim 18, wherein said 
physiological signals include signals indicative of 
hemodynamic information. 

22, A system as in claim 21, wherein said signals 
indicative of hemodynamic information include at least 
cardiac index, blood pressure, pulse oximetry information 
and transcutaneous gas tension. 
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23. A system as in claim 18, further comprising 
obtaining covariate information indicative of a patients 
individual characteristics, and said determining most 
similar patients and said determining information is also 
based on said covariate information. 

24. A system as in claim 21, wherein said 
hemodynamic information includes bio impedance 
information. 

25. A system as in claim 23, wherein said comparing 
and determining most similar patients includes comparing 
temporal patterns within the physiological information. 

26. A system as in claim 23, wherein said 
determining most similar patients comprises 
mathematically determining a state variable indicative of 
the current patient based on at least said hemodynamic 
information and said covariate information. 

27. A method as in claim 23, wherein said 
determining a state variable comprises investigating 
derivatives of the hemodynamic information and integrals 



34 



wo 02/103492 

of the hemodynamic information. 
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28. A method as in claim 18, wherein said comparing 
comprises interpolating among said physiological signals 
to evaluate said physiological signals on a common time 
scale. 

29. A method as in claim 18, wherein said 
determining information comprises modeling said patient 
as a control system with a state parameter representing a 
patient's state, a control input representing a therapy, 
and a nonlinear dynamic system, with process noise being 
explained in terms of covariate information indicative of 
a patients individual characteristics. 

30. A method as in claim 18, wherein said 
determining patients comprises determining quadratic 
distances between a state variable representing a current 
patient and state variables representing other patients 
in said database. 

31. A method as in claim 30, further comprising 
using a weighting matrix to determine said closest 
patients. 
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32 • A method as in claim 31, further comprising 
determining said weighting matrix based on hemodynamic 
parameters and survival versus non survival. 

33. A method, comprising: 

determining hemodynamic parameters of an individual 
patient; 

determining covariate parameters of the individual 
patient representing the patients individual 
characteristics ; 

using said covariate characteristics and said 
hemodynamic parameters to develop a state variable that 
represents information that is continuously variable in 
time; 

comparing said state variable with information 
indicative of a plurality of previously treated patients 
in the database; and 

obtaining survival information about said patients 
in said database. 
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